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Entanglement generation and detection are two 
of the most sought-after goals in the field of 
quantum control. Besides offering a means to 
probe some of the most peculiar and fundamen- 
tal aspects of quantum mechanics, entanglement 
in many-body systems can be used as a tool 
to reduce fluctuations below the standard quan- 
tum limit. For spins, or spin-like systems, such 
a reduction of fluctuations can be realized with 
so-called squeezed states. Here we present a 
scheme for achieving coherent spin squeezing of 
nuclear spin states in semiconductor quantum 
dots. This work represents a major shift from ear- 
lier studies in quantum dots, which have explored 
classical "narrowing" of the nuclear polarization 
distribution through feedback involving stochas- 
tic spin flips. In contrast, we use the nuclear- 
polarization-dependence of the electron spin res- 
onance (ESR) to provide a non-linearity which 
generates a non-trivial, area-preserving, "twist- 
ing" dynamics, analogous to that introduced by 
Kitagawa and Ueda, that squeezes and stretches 
the nuclear spin Wigner distribution without the 
need for nuclear spin flips. 

Recently, squeezing of the collective spin state of many 
atomsi was achieved using atom-light or atom-atom 
interactions^r— . The resulting reduced uncertainty can 
be used to enhance the performance of atomic clocks, 
and to achieve unprecedented precision of measurements 
in atomic ensembles. Similarly, the ability to generate 
squeezed states in a solid-state setting will open new av- 
enues for quantum control of collective degrees of free- 
dom at the nanoscale. In particular, due to the ubiqui- 
tous presence of randomly oriented and evolving nuclear 
spins in nanoscale solid-state devices^"—, future progress 
in spin-based information processing hinges on our abil- 
ity to find ways of precisely controlling the dynamics of 
nuclear spins. 

Recently developed experimental techniques provide a 
powerful toolkit for probing and controlling nanoscale 
groups of nuclear spins. The electron Zeeman splitting 
ofi'ers a sensitive way to detect the collective nuclear 
spin polarization through the Overhauser shift of the 
electronic spin states. This splitting can be measured 
accurately using spectroscopio^iii or time-domain mea- 
surement tcchniques^ i^°'^^i^^ . Furthermore, the hyperfine 
coupling provides a means to control the nuclei through 
their interactions with electron spins. Using such means, 
the Tarucha group has reported nuclear spin polariza- 
tions of up to 40%i^. Using optical pumping, even larger 
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FIG. 1: Nuclear spin squeezing in a quantum dot. a) An elec- 
tron in a quantum dot, with the electron spin S coupled to a 
large group of nuclear spins {In}. Electron spin resonance is 
excited by microwave radiation applied in the presence of an 
external magnetic field, b) Flowchart describing the squeez- 
ing mechanism, c) Schematic depiction of twisting dynamics 
on the Bloch sphere, shown in a rotating frame where the 
mean polarization is stationary. We focus on short to inter- 
mediate times to ^ ^ ^ ii, where the phase space (Wigner) 
distribution is squeezed within a small region of the Bloch 
sphere, with uncertainty decreasing as 1/t beyond the stan- 
dard quantum limit, Eq. (|11I) . At longer times, indicated by 
t2, the distribution extends around the Bloch sphere, at which 
point maximal achievable squeezing is reached. For moderate 
initial nuclear polarizations p = 2I/N, significant squeezing 
can be achieved on the intermediate time scale t < t\ (see 
text). 



polarizations of up to 60% have been demonstrated^^. It 
was proposed that fluctuations of the nuclear polarization 
can be "narrowed" through various schemes exploiting 
electron-nuclear feedback through controlling stochastic 
spin-flips of the nucleii^"— . Recently, such methods have 
been successfully used to achieve narrowing by a factor 
of 5-10 in several systemsiSiiiiHiS^TiSi. These techniques, 
together with the possibility of driving coherent nuclear 
spin rotations using externally applied RF fields^"—, can 
be used to explore a wide range of quantum phenomena 
with nuclear spins. 

Here we describe a spin squeezing mechanism that can 
be employed in a variety of systems such as spins con- 
trolled by microwave radiation in quantum dots deflned 
by electrostatic gates in two-dimensional electron gasesi, 
mesas^, or nanowires^, or optically-controlled spins in 
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self-assembled quantum dots or nanoparticles^rJ^ii^. For 
concreteness, below we provide realistic estimates where 
appropriate using parameters relevant for gate-defined 
quantum dots in GaAs, but we emphasize that the sit- 
uation appears similarly promising for these other types 
of quantum dot systems. 



I. SPIN SQUEEZING IN QUANTUM DOTS 

In a system comprised of many spins {In}, such as a 
quantum dot or an atomic ensemble, the collective to- 
tal spin I = I„ is a quantum mechanical angular- 
momentum variable. Because different vector compo- 
nents of I do not commute, they are subject to the 
Heisenbcrg uncertainty relations 

A/w^>|l(r)l, (1) 

and its cyclic permutations, where A/" = (((5/")^)^/^, 
with 51°' = /" — (/"). The goal of squeezing is to reduce 
fluctuations in one spin component below the "quantum 

limit," A/" = ^\h\{i^)\, set by Eq.© in the situation 
where the fluctuations are distributed equally between 
both orthogonal components. Quantitatively, a state 
is therefore said to be quantum mechanically squeezed 
along an axis z when the squeezing parameteri^ 



is less than one. Because the Heisenberg uncertainty 
principle must remain satisfied, a reduction of fiuctua- 
tions A/^ signified by ^ < 1 must be accompanied by an 
excess of fluctuations A/^ above the quantum limit. 

In contrast to the situation in atomic ensembles where 
classical fluctuations of the collective total spin magni- 
tude can be negligible, however, additional fluctua- 
tions in nuclear spin states in quantum dots are unavoid- 
able due to the abundance of nearly degenerate states 
with energy splittings much smaller than the tempera- 
ture. Therefore, for typical states in quantum dots, the 
uncertainty relation ([T]) is far from being saturated. In a 
typical situation where the equilibrium state oi N ^ 10^ 
microscopic nuclear spin-1/2 moments is completely ran- 
dom, the classical uncertainties A/^ = A/' = ViV^/2 
are identical in magnitude to the quantum uncertainty 
in the maximally polarized (product) state with (/^) = 
hN/2, see Eq.(lT])i^ Below we consider an initial state 
prepared by polarizing nuclear spins to a fraction p of the 
maximal polarization, and then rotating this polarization 
into the equatorial plane of the Bloch sphere such that 
the mean spin points along x, (/'^)o = pNh/2. Because 
the classical fluctuations in the orthogonal components 
of the initial state A/q'^, which are not affected by dy- 
namical nuclear polarization (DNP), are only a factor of 



larger than the quantum mechanical minimum fluctu- 
ations, A/*''^ w only a modest amount of addi- 
tional overhead is required in order to achieve squeezing 
down to the quantum limit. 

We note that, for any mechanism of DNP, the extent 
of polarization p itself will contain fluctuations. Below 
we will check, and conflrm, that such fluctuations do not 
present any significant additional challenges for squeez- 
ing. 

In contrast to recent experiments which achieved 
spin squeezing through atom-light or atom-atom 
interactions-"-, in the mechanism that we outline be- 
low, the nuclear spins arc driven through their interac- 
tion with a single electron spin. As we discuss, the effec- 
tive squeezing Hamiltonian for nuclei is produced after 
motional-averaging of the quickly fiuctuating, strongly 
driven electron spin. While the end results (squeezed col- 
lective spin states) of these approaches are similar, the 
physical mechanisms are quite different. 

II. THE MODEL: SQUEEZING DYNAMICS 

The underlying physics of the feedback mechanism can 
be understood most easily in the regime of fast dephasing 
of the electron spin, which is also a relevant regime for 
current experimental efforts. As depicted in Fig. [T^, we 
consider a single electron in a quantum dot, in contact 
with a large group of nuclear spins, {!„}. The electron 
and nuclear spins are coupled by the hyperfine interac- 
tion IfnF = Sn^"^ ■ where S is the electron spin, 
and each coupling constant An is proportional to the lo- 
cal electron density at the position of nucleus n. The 
electron spin is driven by an applied RF field with fre- 
quency close to the electron spin resonance in the pres- 
ence of an externally applied magnetic field. Because the 
electron spin evolves rapidly on the timescale of nuclear 
spin dynamics, the nuclear spins are subjected to an ef- 
fective hyperfine field (the "Knight field") produced by 
the time-averaged electron spin polarization. Feedback 
results from the dependence of the average electron spin 
polarization on the detuning from the ESR condition, 
which in turn depends on the nuclear polarization, see 
Fig.IB). 

To describe this regime, we model the system with the 
microscopic Hamiltonian (here and below we take ft = 1) 

H = LuzS'-+uJar-+APS' + j{I+S-+rS+)+H,u (3) 

where ujz is the electron Zceman energy in the magnetic 
field, cjQ is the nuclear Larmor frequency, and Hd de- 
scribes the driving of the electron spin and its coupling 
to an environment, which leads to fast dephasing and re- 
laxation. For simplicity, here we consider a single species 
of nuclear spin, and take all hyperfine coupling constants 
to be equal. An = A. The latter condition amounts to 
the assumption that electron density is approximately 
constant inside the dot, and zero outside. In this case. 
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FIG. 2: Time-averaged electron spin polarization S~ , Eq.©, 
and squeezing strength A, Eq.(l7]l, versus RF detuning Jcj from 
the ESR frequency. The average electron spin polarization de- 
pends on P through the dependence of the detuning on the 
Overhauser shift, as indicated by the shaded region. Squeez- 
ing is most efficient for the detuning 5a;* = 7/\/3 where 5*^ is 
most sensitive to P , see text. 
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FIG. 3: Contour plot representation of the Wigner distribu- 
tion of a large collective spin on a locally-flat patch of the 
Bloch sphere, in the rotating frame where wo = 0. The 
mean spin points along x. This picture applies for times 
to < t < t\, see Fig[T]:. Before squeezing, the Wigner dis- 
tribution is isotropic (blue). After squeezing, the Wigner dis- 
tribution, Eg. HlOp . is squeezed along an axis z' , and stretched 
along an orthogonal axis y' (red). Within the mean- field ap- 
proximation, phase space area is preserved. 



the electron spin couples directly to the total nuclear spin 
I = ^„ I„, with the square of the total nuclear spin, 
conserved by the dynamics. The effects of non-uniform 
couplings will be discussed at the end. Due to the large 
mismatch between the electron and nuclear Zeeman en- 
ergies, ujz/ujQ ^ 1, below wc ignore the "flip-flop" terms 
proportional to I~^S~ and I~ S'^ in Eq.([3]). 

We begin by writing the Heisenberg equation of motion 
for the total nuclear spin operator I, dl/dt = i[I,H]: 



K f 



b = 



Wo 



AS' 



(4) 



When electron dephasing is fast on the characteris- 
tic time scale of nuclear spin dynamics, "motional- 
averaging" allows the electron polarization S' in Eq.(|4]) 
to be replaced by an operator-valued scmiclassical mean 
polarization S'{P) which depends on the nuclear po- 
larization P through the Overhauser shift of the ESR 
frequency, c.f. Ref^. In this regime, S'{P) is simply 
determined by rate equations involving the occupation 
probabilities n_ of the up and down electron spin 
states: 



ri+ = W{n^ — n+) + Fin- 



(5) 



Here W = ^n^j/[{6uj- AI')^ is the (/^-dependent) 
ESR transition rate with detuning Suj between the driv- 
ing frequency and ujz, driving strength Q, and electron 
spin dephasing rate 7 = I/T2, while Fi is the electron 
spin relaxation rate. Using this form for W, the steady- 
state solution of Eq. ([5]) gives 

2 {Suj - Apy + 7^ Ti 
where 7 includes the effect of power broadening. 



III. RESULTS: SQUEEZING RATE AND 
EFFICIENCY 



The connection with nuclear spin squeezing is eluci- 
dated by linearizing Eq. (|6]) in AI'^ around the optimal de- 
tuning Jw* = 7/'\/3 where S'^ is most sensitive to nuclear- 
polarization-dependent frequency shifts, see Figl2] Here 
the nuclear Bloch dynamics described by Eq.(|l]) can be 
mapped onto one of the canonical squeezing Hamiltoni- 
ans proposed by Kitagawa and Uedai: absorbing a con- 
stant term into the net nuclear Larmor frequency ujq, we 
obtain 



H KuooP + -X{Pf, \^A 



dP 



1^=0 



(7) 



with P = I{I + 1), I < N/2, conserved by the dynamics. 

The squeezing induced by Hamiltonian ([7]) is most sim- 
ply illustrated by the evolution of a coherent nuclear 
spin state of length / initially oriented along x. For 
large /, scmiclassical analysis shows that the quadratic 
term in Eq. ([7]) induces precession of the total spin vector 
about the z axis with a polarization- dependent Larmor 
frequency r] — dH/dP = + A/^. The resulting evo- 
lution has a "twisting" character, in which an area ele- 
ment on the Bloch sphere is stretched in one direction 
and squeezed in another, see Fig|TJ;. Squeezing becomes 
significant at times t > ts = (|A|/)~^ when the relative 
precession between the upper and lower edges of the cor- 
responding "uncertainty region" becomes comparable to 
the width of the region. Using Eq.ijH]), and assuming that 
the spread of Overhauser shifts AAI' is small compared 
to the width of the resonance 7 to ensure the validity of 
the linearization, we find 



ts 



I6F173 



sVsiA'^-fn^ ■ 



(8) 
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For an order-of-magnitude estimate of the squeezing 
time, we set n = Fi = -17. This choice selects the regime 
of moderately strong electron spin dephasing where the 
resonance is broader than the minimum value 7 = ^Fi. 
In this practically relevant regime, the rate equations ([5]) 
and the motional-averaging approximation can be safely 
applied. Taking the 'intrinsic' width of the resonance to 
be twice larger than the typical Overhauser field fluctu- 
ations, 7 ~ A^/N, wc obtain 



(9) 



t > ts, see Eq.®, the uncertainty A/ of the squeezed 
component decreases as 



tS.min ~ ~f]4~" 



Using a typical value of the hyperfine coupling for GaAs, 
A « 0.1fj,s-\ we obtain ts,min ~ 200 ^is{Vn /I). The 
estimate for ts^min can be improved slightly by optimizing 
expression (jS]) with respect to driving power fi. The fast 
relaxation rate Fi ~ Ay/7^ can be achieved by working 
in a regime of efficient electron spin exchange with the 
reservoirs in the leads. We see that the squeezing time is 
inversely proportional to the initial length of the nuclear 
spin vector, i.e. the degree of nuclear polarization before 
squeezing. 



IV. DISCUSSION: EFFECTS OF 
FLUCTUATIONS 

To investigate the effects of classical fluctuations in 
the initial state, as well as the effect of time-dependent 
electron spin fluctuations around the steady state , we 
now analyze the evolution of the nuclear spin Wigncr 
distribution. For a large initial polarization p, where / = 
pN/2, and for short to intermediate times to < t < ti 
(see FiglT]:), the "uncertainty region" associated with the 
nuclear state is small on the scale of the total spin and we 
can consider evolution in a locally flat patch of the Bloch 
sphere. Here the operators 1^ and approximately obey 
canonical commutation relations, and the initial nuclear 
spin state (polarized along x) is described by an isotropic 
2D Gaussian Wigner distribution with width set by the 
initial transverse fluctuations. A/ = AJq'^. Under the 
polarization-dependent precession generated by Eq.©, 
the Gaussian Wigner distribution evolves to 



[rf + {ly + ixtrf 



2A/2 



(10) 



where without loss of generality we set loq = 0. The ini- 
tial (isotropic) and evolved (squeezed) distributions are 
shown in FigjSl 

The quadratic form in the exponential in Eq. ljlOp is 
diagonalized in a suitably chosen orthonormal basis y', 
z' (see Supplementary Material). As shown in FigOand 
in Eq.(S5) of the Supporting Information, stretching in 
one direction {y') is accompanied by squeezing in the 
perpendicular direction (z'), such that the phase space 
volume of the Wigner distribution is exactly preserved if 
fluctuations of the electron spin are ignored. For times 



A/(t) w A/ 



ts 



(11) 



Squeezing proceeds until long times when the phase space 
distribution begins to extend around the Bloch sphere, 
see FiglT]:. The curvature of the Bloch sphere imposes a 
limit on the maximum achievable squeezing»i 

To derive the squeezing time ts in Eq.®, a coherent 
nuclear spin state with A/css = ■\///2 was used. As 
discussed above, however, when classical uncertainty in 
the nuclear spin state is included, the initial width of 
the Wigner distribution is given by A/ = \/N /2. Given 
that the width A/ of the squeezed component decays 
as 1/t, see Ea. (|TT|) . the effect of the classical transverse 
fluctuations is simply to increase the time required to 
reach a desired level of fluctuations by an order-one factor 

Besides fluctuations in the transverse components of 
the initial polarization, the DNP process used to prepare 
the initial nuclear spin state will also leave behind un- 
certainty in the length / of the net spin (typically with 
a scale much smaller than / itself). Because Eg. pO)) de- 
scribes angular precession with a rate which depends only 
on the z-component of the total spin, however, sections of 
the phase space distribution with constant but vary- 
ing Bloch sphere radii / rigidly preccss without growing. 
Therefore fluctuations in the initial polarization / do not 
pose a significant threat to squeezing. 

In addition to uncertainty in the initial nuclear spin 
state, we must also consider the effect of time-dependent 
fluctuations of the electron spin about its mean-field 
value 5^, Eq. (]§]). The mean- field approximation to 
Eq.Q applies in the motion-averaged limit when the 
electron spin evolves quickly on the time scale of the nu- 
clear spin dynamics, and hence the contribution of time- 
dependent electron spin fiuctuations is small. The resid- 
ual effect of such fluctuations is to add a diffusive com- 
ponent to the nuclear-polarization-dependent precession 
induced by the time-averaged electron spin. As shown 
in the Supplementary Material, the diffusivity n asso- 
ciated with this phase diffusion approximately goes as 
K ~ 1/F, where F '--^ W,Fi is the characteristic rate of 
electron spin dynamics. Thus phase diffusion is indeed 
suppressed by motional-averaging when the electron spin 
evolves quickly on the timescale of nuclear evolution. At 
long times, the competition between coherent twisting 
dynamics, which squeezes fluctuations as and phase 
diffusion, which tends to increase fluctuations as i^/^, 
slows down squeezing to A/ ^ t^^/^, but does not pre- 
vent it. 

The results above are based on a semiclassical mean- 
field treatment of Eq. ^ supplemented by electron-spin- 
fluctuation-driven phase diffusion. This intuitive ap- 
proach is quantitatively supported by a lengthier calcu- 
lation based on the full density matrix of the combined 
electron-nuclear system, to be presented elsewhere^. 
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The more powerful density-matrix approach can also be 
used to study squeezing in the coherent (strong) driving 
regime of electron spin dynamics where the rate equa- 
tions, Eq.(IS]), cannot be applied. 



V. DISCUSSION: EXPERIMENTAL 
FEASIBILITY 

Throughout the discussion, we have worked within the 
approximation of uniform hypcrfinc coupling An = A, see 
Eq. ([3]) , which is widely used for studying electron- nuclear 
coupling in quantum dots. More realistically, hyperfine 
coupling is strong near the center of the dot, where elec- 
tron density is high, and weak at the edges. We note 
that similar variations in coupling occur in atomic en- 
sembles when the size of the atom cloud is comparable to 
the wavelength of light^i22i. Squeezing has been demon- 
strated beautifully in that context, and thus it appears 
that the variation of couplings does not severely impact 
the effect. 

For achievable polarizations of 20%, squeezing sets in 
after ts ^ 2 /is, and fluctuations are suppressed by 
a factor of 10 within approximately 20 fis (neglecting 
phase diffusion). Due to classical fluctuations in the 
initial state, first -^-fold squeezing goes toward reach- 
ing the standard quantum limit, after which quantum 
squeezing proceeds. These timescales are 10 — 100 times 
shorter than the approximately 1 ms-scale nuclear coher- 
ence time recently measured in vertical double quantum 
dots^S. It should thus be possible to squeeze the nuclear 
spin state faster than it decoheres. 

All elements required for achieving and demonstrating 
squeezing have been realized experimentally in nanos- 
tructures at low temperature. Dynamical nuclear polar- 
ization can routinely be produce d^^'^^ and controllably 
rotated using NMR pulses^^"— , and the degree of squeez- 
ing can be ascertained through electron spin dephasing 
mcasurement o^'^'^^i^^'^'^i'^^ . Coherent spin control has 
also been demonstrated in a variety of systemaii^ii^. In 
paricular, we note that in Refi^ electron spin resonance 
was achieved by excitation using microwave magnetic 
fields, with driving amplitudes comparable to the ran- 
dom nuclear field acting on the electron spin, ASI. The 
corresponding transition rates W are of the order of 10 
MHz. In order to reach the motional averaging regime, 
the electron spin relaxation rate, Fi, must be compara- 
ble to the transition rate W, which can be easily accom- 
plished by allowing cotunneling to the electron reservoirs 
next to the dot. 



VI. CONCLUSIONS 

In summary, squeezed states of nuclear spins can be 
produced in quantum dots via feedback provided by their 
hyperfine coupling to an electron spin driven close to res- 
onance. Our squeezing mechanism applies equally well to 



gate-defined dots subjected to microwave driving of the 
electron spin, and to optically controlled spins in self- 
assembled dot ensembles. We have considered various 
physical effects that compete with squeezing, and esti- 
mated the timescale of squeezing. Our estimates indi- 
cate that squeezing is feasible and can be realized with 
current capabilities. Such schemes open the door to un- 
precedented levels of quantum control over collective de- 
grees of freedom in nanoscale systems with mesoscopic 
numbers 10^ to 10^ of nuclear spins. 
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Appendix A: Squeezing of the Wigner Distribution 

In this section, we provide a mathematical description 
of squeezing by analyzing the evolution of a nuclear spin 
state characterized by a Gaussian Wigner distribution. 
As discussed in the main text, for a large spin initially 
oriented in the x direction, and for short times before 
the Wigner distribution extends significantly around the 
Bloch sphere, the Wigner distribution in a locally flat 
patch of Bloch sphere evolves as ft{Iy, Iz) = Ae~^^ '^^ , 
see Eq.(lO), with 

Here A/ = A/q'^ characterizes the transverse fluctua- 
tions in the initial nuclear spin state. 

For times t > 0, the circular Wigner distribution is 
deformed to an ellipse, with major and minor axes deter- 
mined by the quadratic form Q in Eq. (jAl[) . As shown in 
Fig. 3 of the main text, stretching in one direction (?/') is 
accompanied by squeezing in the perpendicular direction 
(z'), such that the phase space volume of the Wigner dis- 
tribution is preserved. The major and minor axes y' and 
z', which lie parallel to the eigenvectors of Q, are rotated 
relative to y and z by an angle 9: 




The angle 9 can be found by extremizing the quantity 
W^wjgw,, w, = ) (A3) 
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with respect to 9. Using the identity [1 — tan^ 0]/2iai\6 
cot 20, we find 



cot 29 = \It/2. 



(A4) 



Note that Eq. (jA4p has two solutions 6*1^2 separated by 
90°, as expected for a symmetric form. 
In the eigenbasis, we write 



ft{P ,r )-^exp 



with 



A4(t) = A/M 1 



1 / ly 



2 \ Al+it) / 2 \ AI-(t) 



(A5) 



(A6) 

In the long time limit A/t > 1, the width AI{t) = A/_ (t) 
of the squeezed component reduces to Eq.(ll). 



by the linear in t twisting/stretching dynamics. There- 
fore for times satisfying t > ^noisc = ^^jtj the coherent 
stretching overwhelms the effect of phase diffusion. 

The efficiency of squeezing in the presence of phase 
diffusion can be estim ated as foll ows. At long times 
t ^ inoiso, the factor + kt/AP describes an increase 
of the width of the Wigner distribution compared to its 
ideal squeezed value AI in Eq.(ll). Combining the t^/^ 
smearing due to phase diffusion with the squeezing, 
we find that the width of the Wigner distribution de- 
creases as at long times: 



A/„ 



AIit)y/l + iit/AP 



XI 



-t 



-1/2 



(B4) 



This expression describes the slowing of squeezing due to 
phase diffusion. 



Appendix C: Calculation of the Phase Diffusion 
Constant 



Appendix B: Phase Diffusion 

The effect of time-dependent fluctuations of electron 
spin polarization about the mean field value can be an- 
alyzed within the rate equation model by introducing a 
time-dependent quantity 



S\t) = S'' + 5S''{t). 



(Bl) 



The fluctuating part can be modeled as delta- 

correlated noise {6S''{t')6S''{t")) oc 5{t'-t"), with an in- 
tensity determined by the rate process, Eq.(5). As shown 
in Supplementary Section [Cj such noise generates phase 
diffusion. 



{S9^{t))^Kt, K^2A- 



{W + Ti)W 
(2W + Ti)^' 



(B2) 



where 69 is the fluctuating part of the Larmor precession 
angle, + iP oc e'^^+''^'. The phase diffusion can be 
accounted for by adding a diffusion term with diffusivity 
k = P K to the equation describing the time evolution of 
the Wigner distribution. 

An important consequence of phase diffusion is non- 
conservation of phase volume, which can be illustrated 
by the evolution of a Gaussian Wigner distribution. Sim- 
ilar to the mean-field case, Eq.(lO), such a distribution 
evolves in time as 



{P-f [P + IXtPf 
2AP ^ 2{AP + kt) 



t > 0. 
. 

Initially, phase diffusion leads to a broadening of 
th e Wigner d istribution, characterized by the factor 
^/1 + kt/AP, which grows like ^1/2 foj. > ^j2^ 

later times, AlXtl > \/ki, the behavior is dominated 



To analyze phase diffusion, we need to calculate the 
generating function for spin fluctuations driven by up- 
down and down-up switching. Denoting the two switch- 
ing rates as W and W', we can obtain the generating 
function for spin fluctuations during the time interval 
< <' < i by approximating a continuous Poisson pro- 
cess by a discrete Markov process with a small time step 
A<t:W-\iW')-\t. We have 



X(A) 



Ra 



N 



1/2 
1/2 



(CI) 



_ , 1 - WA WA 
- 1 w'A 1 - WA 



N 



where W = W + Ti. Taking the limit A 
we obtain an expression 



t 

A' 

0, N 



X(A) = 
M = t 



iX/2 - W 
W 



1/2 
1/2 

W 

~iX/2-W' 



(C2) 
(C3) 



The generating function (|C2p provides a full description 
of the statistics of phase fluctuations, 9t = /J Sz{t)dt, by 
encoding all its cumulants: 



Inx(A) 



00 

E 

k=l 



k\ 



(C4) 



with mi and m2 giving the expectation value {Sz)t and 
the variance {{dt — {9t)Y), respectively. The latter quan- 
tity yields the phase diffusion constant via 7712 = nt. 

Matrix exponential e*^ can be evaluated by writing it 
in terms of Pauli matrices, M = .tq + Xicr;, where 



X3 ==a/2-w_, 

(C5) 
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and we defined W± = {W ± W')/2. We have 



(C6) 



where = x\+x\ + x\ = W^ — X^/i — iXW-. Pluggmg 
this expression for e^^ in Ea. (jC2p . we find 



X(A) = 2e^«* cosh(Xt) + 



smh{Xt) 



-Xl 



(C7) 



an exact expression which is valid both at short times 
and at long times. 

To analyze fluctuations in the steady state, we focus 
on the long times t ^ W~^, {W')~^. In this limit, the 
behavior of x(A) can be understood by replacing coshXt 
and sinhXt by e'^'^, giving 



\nx{X)^iX -W+)t 



X + W+ 



(C8) 



Taylor expanding this expression up to order we find 
the first and second cumulants of phase fluctuations: 



Inx(A) = -iA— — + 



W-t {ixy{wl~wl)t 



/3 



0{X^) (C9) 



Substituting W' = W + Ti,we obtain the time-averaged 
polarization and the phase diffusion constant 



(CIO) 



Crucially, the phase diffusion slows down when the 
switching rates W and W grow, which justifies our mo- 
tional averaging approximation. 
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